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Abstract- In this paper, an optimized model of Shell 
and Tube Parallel Flow Heat Exchanger of water and 
oil type is proposed using C++ programming 
language. Shell and tube heat exchangers are of 
special importance in boilers, oil coolers, condensers, 
pre-heaters etc. They are also used in high pressure 
operations, refrigeration and air conditioning 
industry and process applications. In this paper, a 
number of practical cases of shell and tube heat 
exchangers are taken and the analysis of thermal and 
constructional design of every case is done. The 
optimised design is obtained by minimising the 
pressure drop by maximising the heat exchanger area 
to reduce pumping and running cost. Also 
considering that large sizes lead to increased capital 
cost, heat exchanger area is also optimised to 
overcome this problem. Effect of design parameters 
on pressure drop and heat exchanger area is also 
explained. 


Keywords — Heat Exchanger, C++, Parallel Flow, 
Shell and Tube, Pressure drop 


Abbreviations- 


H.E.= heat exchanger , LMTD= log mean temperature 
difference , P-F= parallel flow, C-F= counter flow, S&T= 
shell and tube, NTU- number of transfer units, LMTD- log 
mean temperature difference, e- effectiveness of P-F HE, 
f= friction factor, Re -reynolds number. 


I. INTRODUCTION 


A. Heat Exchanger:- 


Aheat exchanger is a piece of equipment built for 
efficient heat transfer from one medium to another. The 
media may be separated by a solid wall, so that they 
never mix, or they may be in direct contact. They are 
widely used in space heating, refrigeration, air 
conditioning, power plants, chemical plants, 
petrochemical plants, petroleum refineries, natural gas 
processing, and sewage treatment. 


There are three primary classifications of heat 
exchangers according to their flow arrangement. In 
parallel-flow heat exchangers, the two fluids enter the 
exchanger at the same end, and travel in parallel to one 
another to the other side. In counter-flow heat exchangers 
the fluids enter the exchanger from opposite ends. In a 
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cross-flow heat exchanger, the fluids travel roughly 
perpendicular to one another through the exchanger. 


B. Shell and tube :- 


Two fluids, of different starting temperatures, flow 
through the heat exchanger. One flows through the tubes 
(the tube side) and the other flows outside the tubes but 
inside the shell (the shell side). Heat is transferred from 
one fluid to the other through the tube walls, either from 
tube side to shell side or vice versa. The fluids can be 
either liquids or gases on either the shell or the tube side. 
In order to transfer heat efficiently, a large heat transfer 
area should be used, leading to the use of many tubes. . 
In large power plants with steam-driven turbines, shell- 
and-tube surface condensers are used to condense the 
exhaust steam. 


U-tube heat exchanger 


shell-side 
flue | tube sheet outlet 
baffle shell | plenum 
Y ) 


* 
shell side baffle 
tube bundle - 
with U-tubes shell-side 


fluid out 


Fig. 1.1 SHELL & TUBE HE 
C. Pressure drop in HE- 


HE size can be reduced by forcing the fluid through it at 
higher velocities. But due to higher velocities there will 
be large pressure drops resulting in larger pumping costs. 
The smaller diameter pipe, for a given flow rate, may 
involve less initial capital cost but definitely higher 
pumping costs for the life of the exchanger. 


D. C++ language- 


It is a programming language developed at AT&T Bell 
Laboratories in early 1980s by Bjarne Stroutstrup. ‘C’ 
lacked in simulations and extended the language by 
adding features from Simula 67. The name “C++” was 
coined by Rick Mascitti. It is case sensitive. It supports 
control loops like if , if else, while etc. for conditions. 
The inclusion of header files used for different operations 
are an important part. 
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II. LITERATURE REVIEW In this research, the same case of shell and tube HE of 
water & oil type is taken as in [1] & [2] but for parallel 
A. Previous work- flow and for different design parameters. Multiple 


oo . practical problems are chosen in place of one[1][2] and 
Simplified model of study of thermal analysis of shell analysed using only C++ programming unlike C[1], 
and tube heat exchanger of water & oil type[1][2] is MATLAB[2], ANSYS[1] and FLOEFD[2]. Results 
proposed. It used C code[1] and MATLAB[2] for obtained from thermal and constructional analysis of the 
theoretical calculations of a practical problem of counter multiple problems are used for extracting an optimised 
flow HE. The data obtained is used for designing HE design unlike comparing the results of C[I], 
using Pro-e[1][2]. Also the research include the MATLAB[2] & Pro-e[1][2] with ANSYS[1] and 
comparison of above results with the results obtained FLOEFD{[2]. Also Tc2 and N are a part of given design 
from thermal analysis performed using ANSYS[1] and parameters in the research but in earlier research[1][2], 


FLOEFD[2] softwares. these parameters are calculated . 


B. Scope of work- 


Ii. OBJECTIVES 


TABLE 3.1. UNITS AND SYMBOLS 


PARAMETER SYMBOL SI UNIT 
mass flow rate of hot fluid Mh Kg/s 
hot fluid capacity rate Ch kJ/s °C 
overall heat transfer coefficient between U W/m? °C 
two fluids 
capacity ratio R 
specific heat capacity of hot fluid Cph kJ/kg °C 
mass flow rate of cold fluid Mc Kg/s 
cold fluid capacity rate Cc kJ/kg °C 
specific heat capacity of cold fluid Cpc kJ/kg °C 
rate of heat transfer Q (J/sec) 
length of tube per pass L M 
area of shell Ashell m? 
area of tube Atube m? 
outer diameter of shell Do M 
fluid velocity in tube Vv m/s 
pressure drop P N/m? 
minimum fluid capacity rate Cmin kJ/s °C 
maximum fluid capacity rate Cmax kJ/s °C 
hot fluid temperature at inlet Thl °C 
hot fluid temperature at outlet Th2 °C 
cold fluid temperature at inlet Tcl © 
cold fluid temperature at outlet Tc2 °C 

20 
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diameter of tube 


M 


number of tubes 


In the designing process, the known design parameters 
are- 


e =Thl1= 380°C 
e =©Th2= 360°C 
e =Tel=300°C 
e =Tce2=340°C 


© = Cpe= 4.18 kJ/kg°C 

e Mh=2.5kg/s 

e = Mc=3.5 kg/s 

e = Cph=2.5 kJ/kg°C 

e E=0.00004572 m 

e = U=350 W/m? °C 

e §=(water)=1000 kg/m 

e §—= (oil) =850 kg/m 

A. Design Analysis of multiple cases of HE- 


Various design cases are considered for the thermal and 
constructional analysis. 


TABLE 3.2. VALUES OF d, N AND L 


dl 0.025m 
d2 0.035m 
d3 0.045m 
Nl 20 

N2 25 

N3 30 

Ll 3.5m 
L2 4.5m 
L3 5.5m 


Evaluation of all cases is performed by taking one 
parameter as variable and other two parameters as 
constants from above given values using C++ 
programming. 


B. Study influence of design parameters on As & 
P-The effect of d, N & L on As and P are 
described with the help of Graphical 
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representation performed in MS Excel. The 
graphs are made for the following design cases: 


e When d & N constant- 


1. Lvs As 
2. LvsP 

e When d & L constant- 
1. Nvs As 
2. NvsP 

e When N & L constant- 
1. dvs As 
2. dvsP 


D. Obtain an optimised design- 


The optimum value of design parameters are calculated 
by the program by maximising the HE area. These 
optimum values are then put in the objective function P 
and the other design attributes to obtain the whole 
optimised design. 


IV. METHODOLOGY 
A. Design of HE using thermal and constructional 
formulas- 


The values of known parameters are fed in the following 
thermal and constructional formulas fed in the program. 


e = Thermal formulas- 
Cc=McxCpc 
Ch= Mh xCph 
R= Cmin+Cmax 
Tl= Thl- Tcl 
T2= Th2- Tc2 
LMTD=(T 1 —T2)+ In(T1+T2) 
Q=MhxCphx(Th\-Th2) 
e= (Th1-Th2)+(Th1-Tc1) 
e= (1-{-NTUX(1+R)})/((1+R)) 
NTU= (— In [1-e(1+R)])/(+R) 
e —Constructional formulas- 
Atube =(N1d%2)/4 
Ashell =NadL 
Do = Ashell/nL 
V =(4 Mc)/(1rd’2 QN) 
Re = WVd/(u) 
fl= E/d 
f=[1.1442 log(d/E)] “(-2) 
P= f2QV7L/d 

B. Performing calculations for multiple design cases- 


when Cc>Ch 


opto mo ee op 


where p=0.000596 Ns/m? 


Dame eao op 


The various design cases taken into consideration by 
taking one parameter as variable and other two 
parameters as constants are given below- 
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dl= constant & N1= constant 
d1= 0.025m, N1= 20, L1= 3.5m 
d1= 0.025m, N1= 20, L2= 4.5m 
dl= 0.025m, N1= 20, L3= 5.5m 
d2= constant & N2= constant 
d2= 0.035m, N2= 25, L1l= 3.5m 
d2= 0.035m, N2= 25, L2= 4.5m 
d2= 0.035m, N2=25, L3= 5.5m 
d3=constant & N3=constant 
d3= 0.045m, N3= 30, L1= 3.5m 
d3= 0.045m, N3= 30, L2= 4.5m 
d3= 0.045m, N3= 30, L3= 5.5m 
dl=constant & L1=constant 
d1= 0.025m, N2= 25, L1=3.5m 
d1= 0.025m, N3= 30, L1= 3.5m 
d2=constant & L2=constant 
d2= 0.035m, N1= 20, L2= 4.5m 
d2= 0.035m, N3= 30, L2= 4.5m 
d3=constant & L3=constant 
d3= 0.045m, N1= 20, L3= 5.5m 
d3= 0.045m, N2= 25, L3= 5.5m 
Nl=constant & L1=constant 
d2= 0.035m, N1= 20, L1= 3.5m 
d3 = 0.045m, N1= 20, LI=3.5m 
N2=constant & L2=constant 
dl= 0.025m, N2= 25, L2= 4.5m 
d3= 0.045m, N2= 25, L2= 4.5m 
N3=constant & L3=constant 
dl= 0.025m, N3= 30, L3= 5.5m 
d2= 0.035m, N3= 30, L3= 5.5m 


GpespespespesTpesTpPengTrPensTpPeno espe 


Program of this part is only for calculations. It 
performs the Design analysis of all the cases 
without taking any input. One of the case study is 
shown below: 
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Case study for d1= 0.025m, N3= 30, L3= 5.5m 


C. Study the influence of design parameters on As & P- 


-. 


k. 


So mono sp e 


Thermal design values : 

Ce= McxCpc =3.5x4.21000=14700 kJ/°C 

Ch= Mh xCph= 2.5x2.5x1000=6250 kJ/°C 

Cmin= 6250 kJ/°C 

Cmax= 14700 kJ/°C 

R= Cmin+Cmax =6250+14700=0.4251 

T1= Th1- Tcl =380-300=80 °C 

T2= Th2- Tc2 =360-340=20 °C 

LMTD= (T1 —T2)+ In(T1+T2)=(80-20)= 

In(80+20)=43.2808 

Q=MhxCphx 1000(T/1-Th2)=6250x(380-360 

)=125000 " J/s" 

e= (Th1-Th2)+(Th1-Tc1) 
(380—360)+(380—300)=0.25 

NTU=(—In_ [1-e(1+R)])/U+R)= 

(-In  [1-0.25(1+0.4251))/(1+0.4251)= 0.309 


Constructional design values: 


a. 


Atube = (N7rd’2)/4= (30x7x0.0257)/4 = 
0.01472m? 

Ashell = NardL= 30x1x0.025x5.5= 12.959m? 
Do = Ashell/m = 12.959/(1x5.5) = 0.75m 
V=(4Mc)/(d’2 QN) = 

(4x3.5)/(1x0.025?x 1000x30) = 0.2376m/s 
Re= Vd) = 
(1000x0.2376x0.025)/0.000596 = 9969.43 
fl= E/d = (0.00004572)/0.025 = 0.0018288 
f= [1.14+2 log(d/E)] “(-2) = [1.14+2 
log(0.025/0.00004572)] “(—2) = 0.0228 

p= 
f2QV7L/d=(0.0228x21000*0.2376?x5.5)/0.02 
5 = 567.337N/m? 


TABLE 4.1. GRAPHS AND DESIGN CASES 
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Constant parameters 

Case no. Graph | & 2 Graph 3 & 4 Graph 5 & 6 

1 NILI dill dIN1 
2 N2L2 d2L2 d2N2 
3 N3L3 d3L3 d3N3 
4 NIL2 d1L2 d1N2 
5 N2L1 d2L1 d2N1 
6 N2L3 d2L3 d2N3 
7 N3L2 d3L2 d3N2 
8 NIL3 d1L3 d1N3 
9 N3L1 d3Ll dN1 
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Tube diameter vs Shell area 


= case | 
—— case 2 
====case 3 
== case 4 


=== case 5 


====case 6 


——=case 7 


== case 8 


=== case 9 


CHART 1.d vs As 


Tube diameter vs Pressure drop 


1400 
1200 
= case | 
1000 
—— case 2 
= 800 == case 3 
£ ====_case 4 
@ 
oe 600 == case 5 
====case 6 
400 
==—=case 7 
c a 
aan case 8 
=== case 9 
0 
0.025 0.03 0.035 0.04 0.045 0.05 


d (m) 


CHART 2. d vs P 
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Number of tubes vs Shell area 

25 
20 == case | 
—— case 2 
¢? a See 
= ——case 4 

Dn 

*< 10 = case 5 
=== case 6 
5 === case 7 
== case 8 
===case 9 


CHART 3. N vs As 


Number of tubes vs Pressure drop 

1400 

1200 
== case | 

1000 Fe 

Le — 
= 800 == case 3 
5 = case 4 
a 600 e===ecase 5 
== case 6 

400 
=== case 7 
200 —— case 8 
—SS——— —— 


CHART 4. N vs P 
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Length vs Shell area 


—— case | 
= case 2 
=== case 3 
—— case 4 
— case 5 
=== case 6 


=== case 7 


=== case 8 


= case 9 


CHART 5. Lvs As 


Length vs Pressure drop 


1400 
1200 — case 1 
<~ 1000 —case 2 
= 800 i 
——=case 
= 600 
& 409 ——————o — case 4 
0 ——case 6 
35 4 4.5 5 5.5 6 
——case 7 
L (m) 
CHART 6. L vs P 
D. Obtain an optimised design- e Finding the optimum values of d, N & L 


for which value of P is minimum and 


Optimisation procedure is done as follows: taneimisine HIE! area alse taleite into 


e Selection of minimum value of HE area consideration the capital cost. joe 
from the result data obtained for every Objective — Minimise P= —— 
design case so as to increase its value to the 


e Putting the optimum values of d, N & L in 
other design attributes to obtain optimised 
design. 


maximum. 
e ~=Calculating the values of d, N and L for 
every design case . 
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Vv. RESULTS 
e Simple design of S&T HE is shown in the output screen of C++. 


Temperature of cold fluid(water) at outlet( in degree celsius)=340 
Overall heat transfer coefficient( in Wmetre square degree ius )=350 
Specific heat capacity of hot fluid(oil)( in kJ/kg degree celsius)=2.5 
Specific heat capacity of cold fluid(water)( in kJ/kg degree celsius)=4.2 


Select the number of tubes from the following options 


20 
25 
30 

Enter your choice=20 


Select the diameter of tube( in metres)from the following options 
0.025 
0.035 
0.045 


Enter your choice=0.035 


Select the length of tube Cin metres) from the following options 


3.2 
Enter your choice=< 


Fig. 5.1. INPUT OF DESIGN DATA 


Mass of hot fluid(oil)( in kg/sec)=2.5 
Mass of cold fluid(water)( in kg/sec)=3.5 
Temperature of hot fluid(oil) at inlet( in degree celsius)=380 
Temperature of hot fluid(oil) at outlet( in degree celsius)=360 
Temperature of cold fluid(water) at inlet( in degree celsius)=300 
Temperature of cold fluid(water) at outlet( in degree celsius)=340 
Overall heat transfer coefficient( in W/metre square degree celsius)=350 
Specific heat capacity of hot fluid(oil)( in kJ/kg degree celsius)=2.5 
Specific heat capacity of cold fluid(water)( in kJ/kg degree celsius)=4.2 
Tube diameter (in metres)=0.035 
Tube length Cin metres)=4.5 
Number of tubes=20 
Log mean temperature difference( LMTD)=43 .280849 
Rate of heat transfer( in Joule/sec 
Effectiveness of heat exchanger 
Number of transfer unit 
Area of tube is=0.019233 
Area of Shell is=9.891 
Shell outer diameter is=0.7 
Fluid velocity=0. 181984 
Reynolds Number=10741 .023593 
fr ion factor=0 0956 

sure drop=178 .461281_ 


Fig 5.2. THERMAL & CONSTRUCTIONAL DESIGN OF S&T HE 
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Design table for all cases of d, N & L- 
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TABLE 5.1 DESIGN ANALYSIS 


PARAMETER | dINIL1 | d2N2L2 | d3N3L3 | dIN2L2 | d1N3L3 | dIN1L2 | d1N1L3 

Atube (m’) 0.0098 0.02405 | 0.0477 0.01227 | 0.01472 | 0.0098 0.0098 

Ashell (m?) 5.497 12.37 23.326 8.8357 12.959 7.0685 8.6393 

Do (m) 0.5 0.625 1.35 0.625 0.75 0.5 0.5 

V (m/s) 0.3565 0.1455 0.07335 | 0.2852 | 0.2376 0.3565 0.3565 

Re 14954.1 | 8544.46 | 5538.57 | 11963.3 | 9969.43 | 14954.1 | 14954.1 

fl 0.00182 | 0.0013 0.00101 | 0.00182 | 0.00182 | 0.00182 | 0.00182 

F 0.0228 0.0209 0.01969 | 0.0228 0.0228 0.0228 0.0228 

P (N/m’) 811.356 | 114.0726 | 26.0038 | 669.043 | 567.337 | 1043.17 | 1274.98 

TABLE 5.2 DESIGN ANALYSIS 
PARAMETE | dIN2L1 | dIN3L1 | d2N1L1 | d2N3L3 | d2N2L1 | d2N2L3 | d2N1L2 
R 
Atube (m7?) 0.01227 | 0.01472 | 0.01924 | 0.02886 | 0.02405 | 0.02405 | 0.01924 
Ashell (1m?) 6.8722 8.2466 7.6969 18.1426 | 9.6211 15.1189 | 9.896 
Do (m) 0.625 0.75 0.7 1.05 0.875 0.875 0.7 
V (m/s) 0.2852 0.2376 0.1818 0.1212 0.1455 0.1455 0.1818 
Re 11963.3 | 9969.43 10681.5 | 7121.02 | 8545.23 | 8545.23 10681.5 
fl 0.00182 | 0.00182 | 0.0013 0.0013 0.0013 0.0013 0.0013 
F 0.0228 0.0228 0.0209 0.0209 0.0209 0.0209 0.0209 
P (N/m) 519.26 360.4 138.523 96.488 88.49 139.058 177.626 
TABLE 5.3 DESIGN ANALYSIS 
PARAMETER | d2N3L2 | d3NI1L1 | d3N2L2 | d3N3L1 | d3N3L2 | d3NIL3 | d3N2L3 
Atube (m?) 0.02886 | 0.03180 | 0.03976 | 0.04771 | 0.04771 | 0.03180 | 0.03976 
Ashell (1m?) 14.844 9.896 15.9043 14.844 19.0851 15.5508 19.4386 
Do (m) 1.05 0.9 1.125 1.35 1.35 0.9 1.125 
V (m/s) 0.1212 0.11 0.08802 | 0.07335 | 0.07335 | 0.11 0.08802 
Re 7121.02 | 8307.86 | 6646.29 | 5538.57 | 5538.57 | 8307.86 | 6646.29 
fl 0.00130 | 0.00101 | 0.00101 | 0.00101 | 0.00101 | 0.00101 | 0.00101 
F 0.0209 0.01969 | 0.01969 | 0.01969 | 0.01969 | 0.01969 | 0.01969 
P (N/m’) 79.023 37.064 30.5097 | 16.479 21.1873 | 58.2386 | 37.2896 
27 
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e Optimized design- 
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From the above analysis of design cases, an optimised design is obtained which has optimised value of P and As 
.TABLE 5.4. FINAL DESIGN 


Ce (kJ/°C) 14700 NTU 0.309 
Ch (kJ/°C) 6250 d(m) 0.0262 
Cmin (kJ/°C) | 6250 N 29 
Cmax 14700 L(m) 3.5 
(kJ/°C) 
R 0.4251 Atube (m’) 0.01563 
T1 (°C) 80 Ashell (m?) 8.3544 
T2 (°C) 20 Do (m) 0.7598 
LMTD 43.2808 | V (mis) 0.2238 
Q (J/s) 125000 P (N/m’) 28.3 
& 0.25 

VL. CONCLUSION tube heat exchanger using mat lab and floefd 

software” 

e The optimum value of design parameters [3] Simultaneous optimization models for heat 
d, N & Lis between the maximum and integration by Terrence F. Yee and Ignacio E. 
minimum value hence design is safe. Grossmann 

° The optimum value of other design . [4] Koichi Nakaso, Hiroki Mitani, Jun Fukai 
attributes Ashell, Atube, Do, V & P is “Convection heat transfer in a shell and tube heat 
between the maximum and minimum exchanger using sheet fins for effective utilisation of 
values hence design is safe. eoetey. 

e Pis proportional to L and inversely . 5] Qiu Zhong, Liping Yang, Ye Tao “An optimised 
proportional to d & N? hence decrease in L crossflow plate-fin membrane —based total heat 
and increase in d and N will give a exchanger” 
favourable design. 6] International Journal of Ambient Energy, Volume 

31, Number 4Thermal analysis of counterflow heat 
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